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Abstract 

To achieve accurate results in nanocalorimetry a detailed analysis and understanding of the behavior of the calori- 
metric system is required. There are especially two system-related aspects that should be taken in consideration: the 
properties of the empty cell and the effect of the thermal link between sample and cell. Here we study these two 
aspects for a membrane-based system where heater and thermometer are both in good contact with each other and the 
center of the membrane. Practical, analytical expressions for describing the frequency dependence of heat capacity, 
thermal conductance, and temperature oscillation of the system are formulated and compared with measurements and 
numerical simulations. We finally discuss the experimental conditions for an optimal working frequency, where high 
resolution and good absolute accuracy are combined. 

Keywords: Nanocalorimetry, ac-calorimetry, Membrane, Absolute accuracy, Frequency dependence. 
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1. Introduction 

Thermodynamic measurements with both good abso- 
lute accuracy and high resolution are essential to un- 
derstand fundamental properties of materials. Demand 
for nanocaloric measurements is coming both from the 
wish to study new physics at mesoscopic scales and the 
need to investigate bulk behavior such as anisotropy and 
magnetic field dependence of novel materials that are 
difficult to synthesize as large crystals. A calorimetric 
method particularly suitable for studying sub-|xg sam- 
ples is the ac steady state method (T]-[8]. This method 
usually has high resolution but rather low absolute ac- 
curacy because of the difficulties related to the choice 
of the correct working frequency [9|. A poorly selected 
working frequency results in thermal disconnection of 
the sample, or, for small samples, uncontrolled contri- 
bution of the frequency dependent addenda heat capac- 
ity. A detailed analysis of the behavior of the calori- 
metric system is thus required to know the conditions at 
which the most accurate results are obtained. 

In the primary work by Sullivan and Seidel [ 1 1 heater 
and thermometer were attached directly to the sample, 
which was connected to the thermal bath through a suit- 
able thermal link. Such a design is not feasible in 
membrane-based calorimeters, where both thermometer 
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and heater need to be thin films directly fabricated onto 
a free-standing membrane to maintain the heat capac- 
ity addenda lower than the sample heat capacity. This 
separation of calorimeter cell from sample is not neces- 
sarily a drawback. It simplifies the system analysis and 
increases the experimental reproducibility, since there 
is only one thermal link to the sample that may vary 
from experiment to experiment rather than two or more. 
The most ideal situation is if, furthermore, both heater 
and thermometer are in good thermal contact with each 
other and with the membrane. The thermal diagram of 
such a system is depicted in Fig.fl] Heater, thermometer 
and membrane form a central sample platform, where 
the internal thermal coupling between thermometer and 
heater is assumed ideal. The platform is weakly con- 
nected, through the membrane thermal conductance K e , 
to a thermal bath at temperature 7V The sample is ther- 
mally connected to the platform by a minute amount of 
Apiezon grease or similar. A fundamental requirement 
for accurate measurements is that the thermal link K\ 
between sample and platform is much greater than K e . 
The system of Fig. [T] was analyzed by Velichkov iflOl 
for the case corresponding to a massless thermal link. 
Riou et al. [9| extended the analysis to include the ther- 
mal diffusivity of sample and membrane, following the 
work by Greene et al. ifTTIl . 

In this work we study the system both theoretically 
and experimentally. We show that the system model de- 
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Figure I : Thermal diagram of the studied system. Sample with heat 
capacity C s and temperature T s is coupled through a thermal conduc- 
tance K\ to a platform (Co, To) which, in turn, is connected to the ther- 
mal bath through the supporting membrane. Heater and thermometer 
are thin films lying on top of each other. They compose the platform 
together with the central part of the membrane. The heat capacity Co 
represents an always existing background contribution of the calori- 
metric cell. The membrane outside the platform is treated either as 
a massless thermal conductance K c or as a distributed object with a 
total heat capacity C m . 



scribes our experimental data well, provided that the fre- 
quency dependence of both the membrane addenda and 
the membrane thermal link are taken into consideration 
for small samples. We finally suggest a criterion for the 
working frequency which ensures combined high reso- 
lution and good absolute accuracy. 



2. AC steady state measurements 

In the ac steady state method the temperature of sam- 
ple and calorimetric cell is modulated with a small am- 
plitude, giving rise to a temperature oscillation T ac - 
r ac o sin(a>t - <p) of the thermometer. The power respon- 
sible for this modulation in our case is due to resistive 
heating, and can be expressed as P(t) = Po{l + sinwf), 
where Pq - R^Iq, corresponding to an ac current with 
amplitude Iq y/2 and angular frequency u>' - a>/2 flow- 
ing through a resistor R^. With the use of a lock-in am- 
plifier both the temperature oscillation amplitude r aCj o 
and phase <p are experimentally accessible. They are 
given by flU 



»ac,0 



tan( 



V(wC) 2 + K 2 

coC 

~K 



(1) 



where, for the ideal case of a perfectly connected sample 
and massless membrane, C = Co + C s and K - K e . 
Equation ([TJ can be reshaped to express the unknown C 
and K as a function of the measured parameters (Pq, to, 



C 
K 



wr ac> o 
A) 



sin cf> 
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These expressions form the basics for evaluating ac 
steady state measurements. At low frequency, wt e <k 1, 
where 

T e = (Co + C s )/K e (3) 

is the external time constant, the phase is close to and 
the temperature response is dominated by the thermal 
link. In the opposite limit the phase ideally approaches 
90°. For a>T e > 7 the absolute error arising from us- 
ing the simple expression C = Po/wT- dc ,o is less than 
1%. However, the conditions are rarely ideal at such 
frequencies. To obtain absolute accuracy, a good un- 
derstanding of the system and a carefully selected fre- 
quency are therefore required. 

When measuring small samples, C s < C m , the effect 
of the membrane should be considered. The thermal 
diffusion in the membrane introduces a new frequency 
dependence controlled by the parameter 



C m /K e 



(4) 



For most samples the influence of the thermal link be- 
tween sample and cell is also important. A weak con- 
nection between sample and platform introduces yet an- 
other time scale, described by the internal time constant 



CJKi. 



(5) 



We first study the effect of the membrane, assuming that 
the sample is well connected (r; = 0). We then include 
the effect of a finite thermal link between sample and 
platform. 

3. Effect of membrane 

For small samples, the membrane can no longer be 
approximated as a massless thermal conductance, but 
should be seen as a total heat capacity C m distributed 
over the membrane area. The generated temperature os- 
cillation spreads in the membrane, including the metal 
leads of thermometer and heater, over a distance of the 
order of a frequency dependent thermal length £ t h(oS) = 
y/2D/u> [9 1, where D - K/pc p is the diffusivity, k the 
thermal conductivity, p the density, and c p the specific 
heat. For simplicity, we here model the membrane as a 
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Figure 2: Frequency dependence of r ac o and tan</> for an empty cell 
(C s = 0) at T = 50 K. a) Temperature oscillation amplitude r ac o 
divided by its low frequency value Tfa = Pq/K c . Numerical simu- 
lations are from 1 13 1. b) Phase angle <f> between power and temper- 
ature, expressed as tan0. The model fits are given by Eq. IT) with 
C = Co + C mc ff, K = Kegs, and effective frequency dependences 
from Eq. (6) and 17} . The dashed curves were obtained by using the 
low-frequency limits C = Co + C ra /3 and K = K e . The parameters in 
all cases are Co = 3.18nJ/K, C m = 33.2 nJ/K, and K c = 1.63 U.W/K. 



ID rod. As shown in Appendix A for this model, fol- 
lowing Greene et al. IfTTI . the diffusion practically re- 
sults in an effective frequency dependent addenda heat 
capacity from the membrane, given by 



C„ 



,eff 



C m sinh a — sin a 
a cosh a — cos a ' 



(6) 



with a - y/2(jL>T m . The limits for C m ^s are C m /3 at low 
frequency and C m /a at high frequency. With increasing 
frequency the effective part that is temperature modu- 
lated thus shrinks, decreasing the membrane addenda as 
u) l l 2 . However, not only C = Co + C s + C me ff becomes 
frequency dependent but also K. While the sensed heat 
capacity decreases, the thermal conductance increases 
because of a shorter effective distance of the tempera- 
ture (amplitude) gradients. We show in |Appendix A| 
that it is natural to replace K e by an effective thermal 
conductance 



K e 



eff 



= K P 



a sinh a + sin a 
2 coshff - cos cr 



(7) 



Figure 3: Membrane properties, a) Effective membrane heat capacity. 
b) Effective thermal conductance. The experimental data were evalu- 
ated through Eq. |2). All parameters, experimental data, and assumed 
frequency dependences are the same as in Fig.p] 



The low-frequency limit of K e e (f is again K e , while the 
high-frequency limit is given by aK e /2. 

In Fig. [2] we apply Eq. (J6l) and (J7]i to describe our 
experimental data of a calorimeter cell without sam- 
ple. The calorimeter is made of thin film heater, elec- 
trical insulation and thermometer built on top of each 
other. This stack covers the central 1 10 x 110 \im 2 area 
of a 1 x 1 mm 2 and 150 nm thick Si3N4 membrane, and 
form the platform onto which the sample is placed 0131 . 
The value of K e is easily determined from T^ c at low 
frequency, while C m is obtained from adjusting r m . Fig- 
ure[3]shows the experimentally determined C and K for 
the empty cell, using Eq. d2k, with the corresponding 
model curves obtained by using the same parameter val- 
ues as in Fig. [2] Note that C mie ff is dominating over Co 
at all frequencies. At the highest frequencies, K seems 
to have a tendency to saturate and C decreases slightly 
faster than expected. Whether this is due to an experi- 
mental problem, such as a phase distortion, or if a more 
realistic model would display a different high-frequency 
behavior remains an open question. 



4. Effect of thermal link to sample 

In the normal, experimental case, the effect of a non- 
zero internal time constant limits the upper frequency of 
measurements. The expressions for the stationary tem- 
perature oscillation amplitudes and phase shifts are de- 
rived in Appendix B|for C m = 
and tan 



10' 



0. Surprisingly, r ac , 

are still given by Eq. fJT]) provided that C and 

K are defined by Eq. (|B.10|i and Eq. dBTT). Studying 



these equations, we see that in the low frequency limit 
(coTi <K 1), C — Co + C s and K - K e , i.e., the tempera- 
ture is uniformly distributed in the sample-cell system. 
In the oppsite limit (wr; » 1), the sample is thermally 
disconnected from the platform and just the heat capac- 
ity of the platform is probed. Nevertheless the presence 
of the sample is still sensed by an increase of the appar- 
ent thermal conductance, K — K e + K[. 



From the solution ( B.8 1, we see that the actual sample 
temperature oscillation amplitude decreases faster than 
r ac o when the frequency increases, 



Tac.s _ ^ac,0 V 1 ~~ ■ 



(8) 



where g is defined in Eq. (B.12 1, going from to 1 with 



increasing frequency. The ac signal from a thermometer 
in direct contact with the sample consequently goes to 
zero for <±>t\ » 1, while a thermometer on the cell, as in 
the present case, still gives useful information such as a 
rough estimate of the cell addenda. 

Figure pfl shows T ac o, tan cf> and the transfer function 
in a typical measurement. The sample is a ~ 12.7 \±g 
piece of gold attached to the platform through Apiezon- 
N grease. A clear signature of the decoupling of the 
sample is seen both in T ac o, which no longer goes as 
l/oj, and in tan0, which displays a characteristic de- 
crease during the process. To fit the behavior of Fig. HI 
we introduce the effect of the membrane into the solu- 
tion for a weakly connected sample by replacing Co in 
Eq. (|B.10|i with C + C m>e ff and K e in Eq. dBTT) with 



K e ^ff- We thus have 



C = C + C m , e ff + (1 - g)C s 
K = Xe,eff + gK[ 



(9) 



This is a central result of the paper. The determination 
of the parameters is made easier by first fitting a mea- 
surement of r ac for the empty cell, as shown in Fig. |4^. 
This gives a well-determined value of C m . K e is ob- 
tained from the known power and the low-frequency 
Tacfl, both for the empty cell case and for the case with 
sample. The values of K{, C s , and Co are then ob- 
tained by fitting. A few observations can be pointed 
out. First, the values of K e are slightly different for 
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Figure 4: a) Temperature oscillation amplitude of cell with sample 
and empty cell at T = 227 K. Dashed curves correspond to a well- 
connected sample (g = 0) and no sample (C s = 0). b) larujy for cell 
with sample. Dashed curves correspond to a well-connected sample 
(g = 0) and a disconnected sample (g = 1). c) Transfer function of 
cell with sample. The arrows indicate the frequency where 1% of C s 
has been decoupled (g = 0.01). The fitting parameters of cell with 
sample are: C s = 1.93 pJ/K, K c = 3.09 pW/K, K, = 187pW/K, 
C m = 280 nJ/K, and Co = 42 nJ/K. The parameters for the empty 
cell measurement are X" c = 3.00 pW/K and Co = 25 nJ/K. 



cell with sample and empty cell. This difference is 
within variations between experiments, depending on 
factors such as residual gas conduction and radiation. 
Second, the heat capacity for a Au piece with mass 
m = 12.7 n-g at T = 227 K is about 1.59 |jJ/K from lit- 
erature. The difference between this value and the mea- 
sured heat capacity is coming from the Apiezon grease. 
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Figure 5: a) Heat capacity of cell with sample, b) Thermal conduc- 
tance. Model curves are given by Eq. |9l with the same parameter 
values as in Fig.R] The arrows indicate the frequency where 1% of C s 
has been decoupled (g = 0.01). 



For good absolute accuracy, the heat capacity of the 
grease should thus be measured separately before at- 
taching a sample. Third, there is a difference between 
Co for the empty membrane, and the residual Co when 
the sample is disconnected. This difference can be as- 
cribed to some Apiezon grease in good contact with the 
platform. Fourth, Co of the empty membrane is in our 
case roughly 10% of C m at both discussed temperatures. 
The total membrane addenda can therefore be estimated 
at any frequency provided that the temperature depen- 
dences of C m and K e are known. 

In Fig.[5| we show C and K as obtained from Eq. (|2]i 
and fitted by Eq. d9|, with the same parameters as in 
Fig. El In both Fig.|4]and|5] arrows mark the frequency 
above which more than 1 % of the sample heat capacity 
signal has been lost (g = 0.01 corresponding to wt; m 
0.1). Note that tan0 and K are showing earlier signs of 
the decoupling than r ac o and C. Apparent 1/w behavior 
of r ac o, which corresponds to the plateau of the transfer 
function in Fig. Be, is thus not enough for good absolute 
accuracy. In the middle of the plateau, where experi- 
ments are often performed, the absolute accuracy may 
already be worse than 1%. 



5. Optimizing the working frequency 

In practice, it would be too time consuming to mea- 
sure the frequency dependence of r ac o and tan (f> at ev- 
ery temperature to determine C S (T). It is therefore im- 
portant to establish well-defined measurement condi- 
tions that ensure both good absolute accuracy and high 
resolution. This can be done by adjusting the frequency 
of the measurements, the working frequency « wor k. The 
best w W01 k depends on the requirements and limitations 
of the experiment: absolute accuracy, resolution, al- 
lowed dc offset Tdc, and maximum amplitude of r aCj o. 

The resolution of the heat capacity is obtained 
by differentiating C in Eq. d2|) with respect to 
the temperature oscillation amplitude and phase, 
AC = \dC/dT acfl \ Ar aCj0 + \dCld(j)\ Af For simplicity, 
we assume that the sample is well-connected and that 
the membrane frequency dependence can be neglected. 
We also assume that the phase resolution is given by 
A0 = AT- dc fl/T ac fl, where Ar aC; o is determined by equip- 
ment and setup. Under these conditions 



AC 
~~C 



AT 



ac.O 



1 + 



1 



tanr 



(10) 



At constant r<j c (i.e. Pq), Eq. (10 1 is minimized for 
T RC ,o/Tdc - coscp - 1/ V2 corresponding to <p = 45° 
or tan0 = 1. If we instead maintain a constant r aCj o, the 
resolution can be improved by a factor of 2 by increas- 
ing the frequency further, i.e., by decreasing the excess 
noise factor 1 / tan <p. 

Considering absolute accuracy, lower frequencies are 
generally more accurate as long as Eq. (J2l) is used. Inac- 
curacies (beyond experimental problems etc.) are due to 
a nonzero g in Eq. d5), in turn caused by a finite K{, As 
seen in Fig.Hb, the sample starts to decouple near the lo- 
cal maximum of tan (p. The key parameter that controls 
this maximum is the ratio 



P = KJK e . (11) 

The location of the maximum is <y max Te ~ "VyS, while 
its value is tan0 max » "\^/2. At the maximum (still ne- 
glecting the membrane frequency dependence), the ab- 
solute error is 



1 



1 



P + 2 p 



(12) 



At frequencies a> < oj max , where tan</> » wr e , the abso- 
lute error can be estimated to 



tan<; 



(13) 



With these rules of thumb, consider a few cases. To 
reach tanc/> = 7 (</> max > 82°), as required for using the 
simplified relation C - Po/wT^o with good absolute 
accuracy, a ratio ft > 200 is needed. To have less than 
1 % error at the maximum of tan <p, ft must be greater 
than 100, corresponding to tan</> max > 5 (0 max ^ 79°). 
Finally, to have less than 1% error at tan</> = 1, we only 
need/3 > 10, corresponding to tan^ max > 1.58 (</> ma x ^ 
58°). 



From Eq. (fTOJl and Eq. ( 13 i it is clear that the con- 
ditions are most well-controlled if to is adjusted to 
maintain a constant phase </> wor k. Higher tan</> gives 
less noise, but around tan0 max the accuracy decreases 
quickly and the system will become sensitive to changes 
in ft. Looking at Eq. (13i, it would seem reason- 



able to choose tan as a small fraction of ft. How- 
ever, such a criterion would quickly exceed tan </> max for 
large ft. A better choice is to take tan <p as a fraction 
of tan </> max which goes as yfft. A suitable number is 
tan</> = (2/3)tan</> max , corresponding to 



tan< 



work 



yfft/3. 



(14) 



This condition gives g m 1/(9/3) and thus less than 1% 
error and 100% excess noise as long as ft > 10. If we 
have/3 = 60 as in Fig.|4j we get tan</> =» 2.6 (<f> ~ 69°). 
The expected absolute error at this phase is 0.2% and 
the excess noise is 38%, which is a fairly small number 
considering that a typical achievable resolution AC/C is 
1 in 10 4 to 10 5 . 

Since ft may vary during the measurements, it is im- 
portant to be able to verifying that g remains small. This 
can be done by studying K. At the point where tan </> » 
■\/ft/3 the measured value of K should be ftg » 10% 
higher than K e , which in turn can be obtained from cal- 
ibration measurements. In practice, ft is varying rather 
slowly, and </> wor k can often be maintained at a constant 
value. 

6. Conclusions 

The analysis carried out in this paper illustrates the 
care needed to obtain absolute accuracy in ac steady- 
state calorimeter measurements. Two problems that 
cannot be avoided when studying small samples are the 
frequency-dependent contribution of the sample support 
and the thermal link between sample and support. To 
handle these complications, it is not enough to present 
a model for the temperature oscillation and phase ex- 
pressed as a function of time constants. Rather, explicit 
expressions for the sample heat capacity and external 
thermal link are needed as a function of experimentally 



determinable parameters. Here we provide such expres- 
sions and show that they can be used to overcome the 
experimental obstacles. Based on the analysis, we ar- 
gue that measurements are best performed at a constant 
phase (p. With modern measurement electronics, such a 
condition is both feasible and practical. 
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Appendix A. Membrane frequency dependence 

We approximate the membrane and thin film leads by 
a uniform rod of length L along the x-axis, connecting 
sample at x — 0, at temperature To, with thermal bath at 
x — L, at temperature TV The profile of the temperature 
oscillation can be found from the heat equation 



dT _ d 2 T 
~dt ~ D ~dx 1 ' 



The boundary conditions of Eq. ( A. 1 1 are 



T(0, f) - Tq + r ac , sin cat 
T(L,t) = T h 

where the sample end is subjected to a power 

P = p dc + P(f) = P dc + P sin(wf + 4>). 



(A.l) 



(A.2) 
(A.3) 



(A.4) 



Solving the time-independent problem gives the steady- 
state temperature profile: 



T dc (x) = T + -(T b - T ). 



(A.5) 



The full solution is the sum of the time-independent 
and time-dependent contributions, T(x,t) = T<i c (x) + 
T ac (x, f). The steady-state temperature oscillation along 
the rod is found by means of Laplace transforms: 



T, c (x,t) 



Im I sinh 



f^(x-L) sinhf^zX 



sinh 



I l^-I I 



(A.6) 

This temperature oscillation is connected to the applied 
power through the Fourier law with an additional term 
for the sample heat capacity, 



Pit) 



-LK, 



dr ac (x, t) 



Ax 



+ C SH 



dT, c (x,t) 



.v=() 



df 



(A.7) 



A=0 



where K e is the thermal conductance of the rod and C s +o 
is the sample heat capacity (including central platform). 
Combining Eqs. (|A.2]>, (|A.4|i, (|A.6|l, and (|A.7|i gives 



Inserting the full expressions for P, T s and Ty in 



C s+ 



C m sinh a - sin a 



A. 



a cosh or - cos a wljco 



(A.8) 



and 



^e 



a sinh a + sin a 
2 cosh a - cos a 



COS( 



where the diffusivity D was substituted by L 2 K e /C n 
the parameter a = y/2coT m - 2L/^ t h with r m = C, 



(A.9) 
and 



was introduced. Equation ( A.8 1 and ( A.9 1 take the shape 
of Eq. p) if we define the a-containing terms on the left- 
hand side as the effective heat capacity and thermal link, 
respectively. 

Appendix B. Effect of a weakly connected sample 

In the case of a weakly connected sample, one needs 
to keep track of the temperatures of both sample and 
platform. Below we make the following assumptions: 
(i) the temperature varies so little that the parameters 
are temperature independent, (ii) the thermal links are 
massless, and (iii) the thermal conductances of ther- 
mometer, heater and sample are infinite (as before). 

The thermal response of the platform and sample can 
be written as 

T (t) = T b + T oS + r ac<0 (0 (B.la) 

T s (t) = T h + T oS + T ac , s (f) (B.lb) 

where T ff is the dc offset due to the time-averaged 
power supplied by the heater resistance and any other 
dc power. T ac (t) is the oscillating term that can be fur- 
ther expressed as 

^acoO) = T'acO sin(wf - (p) (B.2a) 

T ac ,s(t) = r ac, s sin(cot - if) (B.2b) 

where <p and ip are the phase shifts which develop be- 
tween power P(t) - Pdc + P{) sin tot and temperature 
T ac (t), due to the finite thermal conductances K e and K[. 
Expressions for the measured variables can be obtained 
with some effort. The thermal equations governing the 
system ensure the conservation of energy for each part 
of the system: 

C0-T7 = P + K e (T b - T ) + ^i(r s - r ), (B.3a) 
df 



dr s 

Cs-r- - Ki(To - T s ). 
at 



(B.3b) 



Eqs. (B.3ai and (B.3bi, we get one time-independent 
equation describing the temperature offset 



T ff - 



r dc 



(B.4) 



and two time-dependent equations, 
coC 



tan(cot - <p) 



+ K e +K { 



Xq -K{X s + Pq sin cot 

(B.5) 

X s —K\Xq 



tan(cot - ip) 
where Xq - T ac sin(o»? - <p) and X s - T acs sin(cot - if) 



Equation (B.5 1 can be expanded through some trigono- 
metric identities before the terms multiplying sin cot and 
cos cot are collected. This results in a system of four 
unknowns: 



(B.6) 



f coCq 


Kz + K, 





-KA 


K e +Ki 


-coCq 


-Ki 








-Ki 


coC s 


Ki 


I Ki 





-Ki 


coCJ 



( X A 




(P6\ 


x 2 







x 3 







A, 




yO) 



where 



fXi\ 

X 2 

x 3 
x, 



fT acfi sin <p\ 
r ac , cos 4> 
r acs sin if 



(B.7) 



\T ac,s COS If) 

The solution to this equation system is 



where 



(Xi) 




( coC 


> 




x 2 

Xi 


= A 


K 

(coC + KconXl - 


-g) 


(B.8) 


J^A, 


\{K - coCcoTi){\ - 


-g» 






A-- 


Po 




(B.9) 




(coC) 2 + K 2 ' 




C = C + (l-g)C s , 




(B.10) 




K = K e +gKi 




(B.ll) 




9 


{con) 2 




m 121 



and 



° 1 + (C0Ti) 2 ' 

The solution given in this form can be easily reshaped 
into Eq. (fTh for the experimentally measured values J a c,o 
and <p if C and K are defined as above. 
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